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We investigate the impact on the classical dynamics of dark matter particles and dark energy 
of a non-minimal coupling in the dark sector, assuming that the mass of the dark matter particles 
is coupled to a dark energy scalar field. We show that standard results can only be recovered if 
the space-time variation of the dark energy scalar field is sufficiently smooth on the characteristic 
length scale of the dark matter particles, and we determine the associated constraint dependent on 
both the mass and radius of the dark matter particles and the coupling to the dark energy scalar 
field. We further show, using field theory numerical simulations, that a violation of such constraint 
results in a microscopic feedback effect strongly affecting the dynamics of dark matter particles, 
with a potential impact on structure formation and on the space-time evolution of the dark energy 
equation of state. 


I. INTRODUCTION 


Cosmological observations provide overwhelming evi¬ 
dence that our Universe is currently undergoing an in¬ 
flationary phase m- In the context of Einstein’s Gen¬ 
eral Relativity, this accelerated expansion can only be 
explained if the Universe is presently dominated by an 
exotic Dark Energy (DE) form, violating the strong en¬ 
ergy condition [zhu]- While DE explains the observed 
dynamics of the Universe on cosmological scales, a non- 
relativistic Dark Matter (DM) component is required in 
order to account for the observed dynamics of cosmolog¬ 
ical perturbations over a wide range of scales. Together 
DM and DE seem account for about 95% of the total 
energy density of the Universe at the present time. How¬ 
ever, our knowledge of DM and DE is indirect, relying 
only on their contribution to the gravitational field. Un¬ 
veiling the nature of DM and DE is therefore one of the 
most ambitious challenges of fundamental physics. 

The possibility of a non-minimal coupling in the dark 
sector oma is an exciting topic of current research. 
Such a coupling could affect both the background evo¬ 
lution of the Universe as well as on the growth of cos¬ 
mological perturbations |16fI20| . with a potential im¬ 
pact on the redshift dependence of the apparent magni¬ 
tude of type la supernovae, the Cosmic Microwave Back¬ 
ground (CMB) anisotropies, baryonic acoustic oscilla¬ 
tions, redshift-space distortions and weak gravitational 
lensing (see, e.g., m and references therein). A non- 
minimal coupling between DM and DE may naturally 
lead to a fractional matter abundance at the present time 
significantly different from the one obtained using the 
information contained in the CMB temperature power 


spectrum assuming no interaction between DM and DE 
[22] , which can help resolving the tension between CMB 
and local constraints on the value of the Hubble param¬ 
eter. In fact preliminary indications of a late-time non- 
minimal interaction between DM and DE have been re¬ 
cently reported in (see also m- The coupling 

of DE with other fields has also been investigated in the 
context of varying fundamental constants (see, e.g., EH) 
and growing neutrino models (see, e.g., [2BJGS3)- 


In previous studies, the coupling between DM and DE 
energy has been considered at a macroscopic level. These 
studies (see, e.g., EH 1301) implicitly assume the space- 
time variation of the DE field to be sufficiently smooth 
on the characteristic length scale of the DM particles, 
thus neglecting any non-linear feedback at the micro¬ 
scopic level. In this paper we intend to bridge this gap 
by studying the impact of a non-minimal coupling in the 
dark sector on the microscopic dynamics of DM particles 
and DE. In Sec. [B]we start by presenting a simple field 
theory model, where the mass and size of the DM parti¬ 
cles is assumed to be a function of a DE scalar field, and 
we use it to derive the standard results for the dynamics 
of DM particles non-minimally coupled to the DE field as¬ 
suming a negligible microscopic feedback. In Sec. Ill we 


perform field theory numerical simulations in 1 + 1 dimen¬ 
sions and quantify the effect of the microscopic feedback 
neglected in the previous section. We discuss the poten¬ 
tial cosmological implications of our results and present 
the main conclusions of this work in Sec. hyi 
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Throughout this paper we use units such that c = 1, 
where c is the value of the speed of light in vacuum, and 
we adopt the metric signature (—, +, +, +). 
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II. INTERACTING DARK ENERGY MODEL with 


In this paper we study the microscopic dynamics of 
non-minimally interacting DM and DE. Although the 
main results derived in the paper are generic, for sim¬ 
plicity we shall assume DM and DE to be described by 
two coupled scalar fields, <j> and p, living in a 1 + 1 di¬ 
mensional space-time. We shall consider a class of models 
described by the action 

s = jd 2 xV^gC 1 (l) 

where the Lagrangian C is given by 
£ = ~d l j(t>d ll <l>-'^d v <pd , '(p-V(<l>,<p). (2) 


Here, g = det(g^) and g^ v are the components of the 
metric tensor. For concreteness, we shall assume the fol¬ 
lowing form of the potential 

V(4>,v) = v (t,v) + W(v), (3) 

= ^V-i 2 ) 3 , (4) 


where the scalar fields 0 and p play essentially a DM and 
DE role, respectively, with the non-minimal interaction 
between them being realized through the A = X(p) term. 
The equations of motion for the scalar fields 


□0 = 



(5) 


\3<p 


dW 

dp 


m, 


( 6 ) 


may be obtained by minimizing the action with respect 
to variations of 0 and p. Here, □ = V M V M , V M represent 
the covariant derivates and 


_ 1 din A 

P ~ ~ 2 dp ' 

The energy-momentum tensor 

T^ = _JL = 2-^- + g^C 

y~9 oQlli' 

is given by 


(7) 

( 8 ) 


Tllv = t dm + t de > (9) 

with 

a* iV 

T% m = ^0cU0 - y — (3 a 03“0 + 217(0, p)) , (10) 

a^ 

T%e = d^pd v p - 9 — (d a pd a p + 2 W{p)) . ( 11 ) 

The energy-momentum tensors of the DM and DE com¬ 
ponents are not individually conserved since 

= Q* , ( 12 ) 

= -Q", (13) 


= 2/317(0, p)d li p. 


(14) 


However, the total energy-momentum tensor is covari- 
antly conserved (V^T^^ = 0). Further ahead, we shall 
demonstrate that Eq. (12), ignoring microscopic feed¬ 
back, leads to the following equation for the dynamics of 
the DM particles 


dn a 

— = (V/3 p a ) u 13 = m<i > l3d a p . (15) 

ar 

where p a = m^u 01 , and r and u a are, respectively, 
the mass, the proper time and the coordinates of the 
N + 1-velocity of the DM particle (here N represents the 
number of spatial dimensions). 


A. Dynamics in 1 + 1 dimensions 


In this paper, we investigate the dynamics of DM 
particles non-minimally coupled to DE, focusing on mi¬ 
croscopic feedback with a characteristic timescale much 
smaller than the cosmological timescale. Therefore, in 
the following, we shall neglect the expansion of the Uni¬ 
verse. In a 1 + 1 dimensional Minkowski space-time the 
line element can be written as ds 2 = — dt 2 + dz 2 so that 
the equations of motion for the scalar fields 0 and p are 
given respectively by 


0-0" 


p-p 


dU 
30 ’ 
dW 
dp 


+ 2 pu. 


(16) 

(17) 


where a dot denotes a derivative with respect to the phys¬ 
ical time t and a prime represents a derivative with re¬ 
spect to the space coordinate 2 . 

The components of the energy-momentum tensor of 
the DM and DE components can now be written as 


^7-i00 

1 DM 

rjiOz 
1 DM 


rpZZ 

1 DM 


rpOO 
1 DE 

rpOz 

1 DE 


rpZZ 

1 DE 
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+ £7(0, p ), 


A2 ±12 



— ~U{(j>,p), 

p' 2 

\ + W(p ), 

5 

p' 2 

\ - W(p ), 


(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


Let us start 
assume that 0 
(16) becomes 


by considering a static DM particle, and 
= 0(z) and A = const. In this case Eq. 


dU 
d0 ’ 


( 24 ) 
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and it can be integrated to give 



(25) 


assuming that \<j)\ —> p for z —> 
that in the static case the total 
to 2U. If the particle is located 
the following solution 


±oo. Eq. (25) implies 
energy density is equal 
at z = 0, Eq. (24) has 


cA = ±p tanh 



(26) 


The equations of motion (16 1 and © are solved using 
a 2nd order Runge-Kutta time-integration method com¬ 
bined with a 2 nd order finite difference scheme for the 
spatial derivatives. In all the simulations we assume that 


ft — /^onO( + Pon ) /^onO(+ Poff) > 


(32) 


where 0 is the Heaviside function. This implies that /3 
is switched on to /? = (3 on for p = p on and switched off 
for p = tpoff- Using Eqs. (|29| and (32) one finds that 


with 

R = A - 1 / 2 ? 7 -1 . (27) 


The rest mass of the particle is given by 

nOO 

T 00 dz = 2 


TO 0 = 


8 

Udz — - U rnax R — 


2 y ^2 


A 


which implies that 
p _ 1 din A din m$ 


2 dip 


dip 


(28) 


(29) 


Here U max = U(4> = 0) = Ar? 4 /4. 

The evolution of the total energy E = and linear 
momentum p = associated with the DM particles 

may be calculated using Eqs. (10) and (12) 


(i § ~ f d ° T DM dz - 

= 2 J /3U(<j>, <p)d°(pdz ~ — fi m <tP , (30) 

F =f t = ~J d ° T D M dz = 

= 2 J /3U(<j>, <p)d z ipdz ~ , (31) 


Hon 


1 ln ( m 4Poff) \ 

tpoff tp on \ TTi(j){fp on) ) 


(33) 


Our simulations run in a one-dimensional grid of size 
L = 1 x 10 4 with periodic boundary conditions. The 
initial conditions for the scalar fields (</> and p) and their 
time derivatives are given by 




Pi 

Pi 


p ^tanh 
tanh 


7i(z- £/4) 
\/ 2 Ri 

li{z - 3L/4) A 


V2Ri 

±v i 'Yi*f2U((f)i,ipi ). 




0 , 

4 x 10" 3 , 



(34) 

(35) 

(36) 

(37) 


where v t = 1 x 10~ 3 , Ri = 20, the subscript ‘i’ refers to 
the initial time U = 0, and the + sign in Eq. ( |35| ) applies 
for z > L/ 2, while the — sign applies for z < L/2. The 
initial conditions for (f> and <f> describe a DM particle and 
an antiparticle at a distance L/2^> R from one another, 
moving both with velocity v = v l . The final time t = tf = 
5 x 10 3 is such that the particle and the anti-particle are 
effectively isolated from each other during the course of 
the simulation (the fifth force between them, mediated 
by the DE scalar field, would only be felt for t > tf). 

In our simulations we take 


where v is the DM particle velocity and 7 = (1 — w 2 ) -1 / 2 . 
Here, we assume that ip and p' are smooth on the charac¬ 
teristic length scale of the particle, which is not the case 
in general. Eq. (31) implies that, in the absence of gradi¬ 
ents of the DE field p, the DM particles do conserve their 
linear momentum. The gradients of the field are respon¬ 
sible for a fifth force mediated by the DE field p which, 
in the Newtonian limit, is / 3 2 /(47 tG) times stronger than 
gravity. Taking into account that p° = E, p z = p , and 
dt/dr = 7 , it is simple to show that the covariant form 
of Eqs. (30) and (31) is given by Eq. (15). 


III. NUMERICAL SIMULATIONS 

In this section we describe the results of field theory 
numerical simulations in 1 + 1 dimensions of the dynam¬ 
ics of DM particles interacting non-minimally with DE. 


£ = m^ipof^/m^ipon) , 


(38) 


to be either £ = 1.2 or 0 . 8 , corresponding to a model 
where the mass of the DM particles increases or decreases 
by 20%, respectively. Our units of velocity and energy 
are such that c = 1 and p = 1 (except for these conditions 
the choice of units is arbitrary). For simplicity we shall 
also assume that W = const. 

Let us define the parameter a by 


('^°// ^on)p Ri 

a = -— „ —-— X - 

2 Pi A t v 


(39) 


where A t v is a characteristic time associated to the vari¬ 
ation of p defined as 


A t v = (Poff 

Pi 


(40) 
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FIG. 1: The top and middle panels show, respectively, the 
evolution of f EM and /dm (solid lines) and of f EE and f EE 
(dashed lines) considering £ = 1.2 for a — 2 (black and red 
lines) and a = 0.2 (blue and magenta lines). Note that the 
values of f° and f z , represented by dotted lines in the top and 
middle panels, respectively, are always very close to unity as 
required by energy conservation. The lower panel shows the 
evolution of r u and r G for a — 2 (black and red lines, respec¬ 
tively) and a = 0.2 (blue and magenta lines, respectively). 


In our simulations we set ip on = 1 x 10 2 </ji and p 0 ff is de¬ 
termined by the value of a , which we choose to be either 
0.2 or 2. The parameter a represents the ratio between 
the average variation of the characteristic microscopic en¬ 
ergy density of the particle on a timescale equal to Rt, 


AU max Ri/ (Af v ), 

and the initial kinetic energy density associated with the 
DE, ifii/2. For a > 1 microscopic feedback effects associ¬ 
ated with the evolution of the mass of the DM particles 
are expected to be large. On the other hand, for a -C 1 


FIG. 2: Analogous to Fig. |T] for £ = 0.8. 


these effects are expected to be small and the standard 
results for the dynamics of DM particles non-minimally 
coupled to the DE field should hold. 

The following functions 


/dm = 

/dm = 
/° = 
r G = 


_ I Tdm d z 


ro — f 1 d°e dz 


foE = 


fi>E = 


fT t °2 0 dz ’ 
jT°R M dz 
f 1't=o dz ’ 

/dm + Ide 
/ </' 2 dz 

2 / ^DM|t=0 dz ’ 


fT t °° 0 dz’ 

I Tde dz 
f Tt=o dz ’ 
f z = /dm + Ide 5 
r u = JUdz 


r T oo ’ 

J "*■ DM\t =0 


(41) 

(42) 

(43) 

(44) 


were computed with our simulations for £ = 1.2 (Fig. 
[l]) and £ = 0.8 (Fig. [2]). Here we have disregarded the 
contribution of the DE potential W (ip) in the calculation 
of and T 00 . 

Fig. [I] (top panel) shows the evolution of fp M (solid 
lines) and fp E (dashed lines) for a = 2 (black and red 
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lines) and a = 0.2 (blue and magenta lines). It shows the 
increase of f EM by 20%, associated to the increase of the 
mass of the DM particles, as well as the corresponding 
decrease of f EE by the same amount. The evolution is 
faster in the a = 2 case than for a = 0.2. The value of 
/° (dotted line) is always very close to unity, both in the 
a = 2 and in the a = 0.2 cases, as required by energy 
conservation. 

Fig. [l] (middle panel) shows the evolution of f EM 
(solid lines) and f EE (dashed lines) for a = 2 (black 
and red lines) and a = 0.2 (blue and magenta lines). It 
shows that the decrease of f EM , associated to the de¬ 
crease of the linear momentum of the DM particles, is 
exactly compensated by the increase of f EE , so that f z 
(dotted line) is always very close to unity, as required 
by linear momentum conservation. For a = 2 the lin¬ 
ear momentum of the DM particles is changed dramat¬ 
ically while for a = 0.2 it is nearly conserved. This is 
not surprising, since microscopic feedback should have a 
negligible impact on the dynamics of the DM particles 
for a 1, only becoming significant for a > 1. 

Fig. 0 (lower panel) shows the evolution of r u and r G 
for a = 2 (black and red lines, respectively) and a = 0.2 
(blue and magenta lines, respectively). It shows that 
r u and r G have always a similar magnitude (note that 
r u = r G in the static case) and present long lived oscil¬ 
lations (more pronounced for a = 2 than for a = 0.2). 
Notice that r G and r u have a phase lag of 7r and that 
r u + r G evolves from approximately 1 to 1.2 as a di¬ 
rect consequence of the 20% decrease in the mass of the 
DM particles. The DM energy associated with the time 
derivative of <j) increases with a but it is very small both 
for a = 0.2 and a = 2. 

Fig. 0 is analogous to Fig. |T] but now for £ = 0.8, 
meaning that the mass of the DM particles decreases by 
20%. One of the most significant differences between 
Figs. [I]and[2]is related to the very different timescales for 
the energy-momentum transfer between the DM particles 
and the DE held for a = 2, which happens even though 
the value of A t v is the same for £ = 1.2 and £ = 0.8. This 
has to do with the fact that a rapid increase (decrease) of 
the mass of the DM particles is also associated to a rapid 
decrease (increase) of the kinetic energy of the DE scalar 
held, thus leading to a switch off time signihcantly larger 
(smaller) than A t v for £ = 1.2 (£ = 0.8) in the a = 2 
case. Hence, for a > 1 the microscopic and macroscopic 
values of (3 can be in general be very different in models 
where /3 = (3(ip). 


sidered only at a macroscopic level, and it should eventu¬ 
ally be validated at a quantum level. We have shown us¬ 
ing analytical arguments and numerical simulations that 
although the standard results for the dynamics of the DM 
particles can be recovered if the space-time variation of 
the DE scalar held is sufficiently smooth on the charac¬ 
teristic length scale of the DM particles, in general the 
DE coupling introduces non-trivial backreaction effects. 

We have found that microscopic feedback becomes im¬ 
portant when the variation of the characteristic micro¬ 
scopic energy density of the particle on a timescale equal 
to R is of the order of kinetic energy density associated 
with the DE. In three spatial dimensions this is equiva¬ 
lent to the condition 

+ 1 ). ( 45 ) 

where p v and p v = w v p v are the DE energy density and 
pressure. Hence, microscopic feedback should be taken 
into account for DM masses 


> m* 


P<p{wtp + 1 )R 2 

pip 


(46) 


As a reference, a DM particle with a radius and mass 
similar to that of the proton ( m p ~ 10 ~ 2 ‘ kg and r p ~ 
10~ 15 m) microscopic feedback would be important today 
even for a value of (3q as small as /?o/>o ~ 10~‘ l H o (here a 
‘ 0 ’ refers to the present time and we have assumed that 
u^o — 1| < 0.1). Note also that, m* —► 0 for R —>• 0, 
irrespectively of the value of (3. 

We have demonstrated that for > to* a signifi¬ 
cant transfer of linear momentum between the DM par¬ 
ticles and the DE scalar field is associated to microscopic 
feedback, which could have an important impact on the 
growth of cosmic structures. On the other hand, for 
TO^ <C to* the linear momentum transfer due to micro¬ 
scopic feedback is negligible. We have also found that 
for ui^ > to * the microscopic energy and linear momen¬ 
tum DE density around the DM particles can in general 
be very different from the average ones. This implies 
that microscopic feedback may play a relevant role on the 
space-time evolution of the DE equation of state, and can 
be associated with very different microscopic and macro¬ 
scopic values of (3 in models where (3 = /3(<p). The re¬ 
quirement that microscopic feedback effects do not have 
a dramatic impact on the dynamics of both DM and DE 
may turn out to be one of the strongest constraints on the 
nature of DM particles in coupled dark energy scenarios. 


IV. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the classical micro¬ 
scopic dynamics of DM particles non-minimally coupled 
to a DE held, using a simple model which is expected to 
capture the essential features of generic scenarios where 
the mass of the DM particles is coupled to the DE scalar 
held. This analysis complements previous studies where 
the coupling between DM and DE energy has been con¬ 
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